UCRL-JRNL-215137-DRAFT 



Exploring Lifetime Effects in Femtoscopy 

D.A. Brown, R. Soltz. and J. Ncwby 
Lawrence Livermore National Laboratory, Livermore California 94-551 

A. Kisiel 

Faculty of Physics, Warsaw University of Technology, ul. Koszykowa 75 00-661 Warsaw, Poland 

(Dated: February 1, 2008) 

We investigate the role of lifetime effects from resonances and emission duration tails in femtoscopy 
at RHIC in two Blast- Wave models. We find the non-Gaussian components compare well with 
published source imaged data, but the value of Rout obtained from Gaussian fits is not insensitive 
to the non-Gaussian contributions when realistic acceptance cuts are applied to models. 
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The experiments at the Relativistic Heavy Ion Collider 
(RHIC) have produced a series of remarkable results in- 
cluding the discovery of jet-quenching through the sup- 
pression of high particles [1] and the observation of hy- 
drodynamic flow [2] . The complete set of measurements 
point to the existence of a dense partonic medium [3- 
6] that evolves hydrodynamically from the time of ther- 
malization until freeze-out. Several hydrodynamic mod- 
els have succeeded in reproducing measured spectra and 
flow [7-9], but with the exception of parameterized fits 
to the data [10], they have not reproduced the space-time 
measurements for the particle emission region [9, 11-15]. 
While some the models match the systematic trends of 
the longitudinal Gaussian extent, Riong, they produce an 
extent which is too large in the transverse component of 
the particle pair emission direction, i?o„(, too small in the 
orthogonal transverse direction, Rside, or both. The ra- 
tio of these two measurement, Rout /Rside, has been sug- 
gested as an indication of the duration of the freeze-out 
emission stage [16], and its near unitary value for a wide 
range of colliding systems, geometry, and center-of-mass 
energy in light of the disagreement with hydrodynamic 
models has been called the "HBT Puzzle" . 

The traditional spatial-temporal analysis of the parti- 
cle emission region in a high energy collision pioneered 
by Goldhaber, Goldhaber, Lee, and Pais [17] assumes a 
parameterized source shape, usually Gaussian, which is 
fit to a measured enhancement in the two-particle mo- 
mentum distribution at low relative momentum. The 
physical basis for this technique is analagous to intesity 
interferemetry techniques developed by Hanbury Brown 
and Twiss (HBT) [18] for which the technique is often 
named. It has been applied mostly to pairs of identi- 
cal pions emitted in both lepton and hadron collisions, 
but the greatest interest has been generated in the heavy 
ion collisions, where space-time measurements hold the 
promise of providing constraints on the nuclear matter 
equation of state. The technique has since advanced in 
many ways, incorporating a variety of systematic depen- 
dencies on pair momentum, collision geometry, reaction 
plane and a decomposition of the source emission region 



into three cartesian radii, i?side, -Rout, -Riong and cross- 
terms. However, the need to parameterize the source 
distribution by its "lengths of homogeneity" [23, 24] has 
remained a fundamental limitation (see [19] for a recent 
review and references therein). This limitation is un- 
derscored by the application of three-pion correlations 
to rule out coherent pion production as an explanation 
for the non-unitary value of the correlation strength. A, 
which provides evidence for a complex source shape con- 
sisting of core and extended halo [20-22]. 

The tools to resolve non-Gaussian components have re- 
cently been provided by a source imaging technique [25- 
27] and the first evidence of non-Gaussian tails have been 
observed by the PHENIX experiment in 200 GeV cen- 
tereal Au-I- Au Colhsions [28] , thereby demonstrating the 
ability to image non-Gaussian sources and resolve contri- 
butions from an extended lifetime and/or resonance halo. 
Understanding the role of these two contributions to an 
extended source distribution has important consequences 
for the "HBT Puzzle" and the limitations of traditional 
techniques. In this paper we will explore these limitations 
and their implications using two models to decouple the 
relative contributions from lifetime effects and resonance 
decays. 

Our first model uses a blast-wave fiow profile for the 
source distribution [32] , but substitutes exponential time 
emission and Gaussian longitudinal profiles for a Gaus- 
sian in proper time and an infinite Bjorken tube. This 
model also includes a resonance contribution from uj feed- 
down. 

The building blocks for our two-particle source arc the 
normalized particle emission rates for the core and halo: 

Dir, t, p) - /i?corc(r, t, p) + (1 - ./) Aiaio(r, t, p) (1) 

The core component (fraction = /), consists of a Gaus- 
sian spatial part, an exponential time profile, and a mo- 
mentum dependence arising from hydro-like Boltzman 
factor: 

(2) 
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Here, /, T, R^, Ry, i?^, and r/o are all adjustable and 
are specified in the lab frame. For this study wc set 
T = 165 MeV, / = 0.5, = Ry = Rt ^ Rz = 4, fm, 
and Tfo = 10 fm/c. The flow profile is given by 

(^) = (cosh 7] cosh p, fx sinh p, sinh 77 cosh p) (3) 

where ry = i In (f^f ) and p = po'^rZ-Rmax with po = 0.6. 
These parameters were chosen arbitrarily since this paper 
is a generic study of the time profile and lj contributions 
to the source distribution. For comparison, the Blast- 
Wave fits of Retiere and Lisa [32] for central collisions 
converged for R^ « Ry w 13 fm (equivalent to a 2D rms 
Gaussian radius of 6.5 fm), T w 110 MeV, and po ~ 0.9. 

Following earlier work on the contribution of resonance 
decays to the pion source distribution [35, 36], we as- 
sume that the halo is dominated by the decay of the 
u) resonance {t^ = 23 fm/c), with a fractional contri- 
bution of 50% to the pion distribution in the region of 
0.2 < kx < 0.36 GeV. Other potential candidate reso- 
nances have decay times that are much too short, e.g. 
the p with lifetime Tp = 1.3 fm/c, or too long, e.g. the 
77' with lifetime r^' = 975 fm/c or have charged pionic 
decay modes with small branching fractions. The w's are 
emitted from the same core, but we allow them to prop- 
agate classically for some distance before decaying into 
pions: 

£'haio(r, t, p) « / dM P(p.., p)e"^*/"- 

^ p (4) 

xDeorc(r- -^At,t- At, p^). 

We include the dominant three-body reaction lo — > 
7r+7r^7r° with a branching fraction of 88.8% and define 
the probability, P, for finding a tt with momentum p from 
the decay of the lo with momentum Pi^, using standard 
three-body decay kinematics. We neglect the only other 
charged pion decay mode of the u which is w — > ■k'^'k^ 
with a branching fraction of 2.2%. This mode would only 
have a negligible impact on the source distribution. 

This source function for the probability to emit a pair 
with a separation of r' in the pair CM frame is given 
by [25] 

5p(r') = j dr'^ J d^RD,{R+r/2,P/2)D2{R-r/2,P/2). 

(5) 

We construct the source function from this single parti- 
cle source by Monte-Carlo integrating our emission func- 
tion in Eq. (5). We work in Bertsch-Pratt coordinates 
[30], so that the time integral in Eq. (5) serves to move 
the time effects into the outward and longitudinal di- 
rections. To compare to the ID, angle-averaged source 
image for central Au-Au collisions recently measured by 
PHENIX [28], we keep only pairs where both pions have 
a pseudo-rapidity \r]\ < 0.35 and transverse momentum 



0.2 < kr < 0.36 GeV. We also applied the STAR ac- 
ceptance cuts in pseudo-rapidity, [ry] < 1.0, but found no 
qualitative difference. 

In Fig. 1, we plot our "basehne" model source func- 
tion as 1-fm slices along the various directions in the 
Bertsch-Pratt coordinate system. The top three pan- 
els show slices for Rside, Rout, and Riong, respectively, 
and the bottom panel shows the full angle-averged source 
{£m = 00 term in a spherical harmonic expansion of the 
source). All distributions are normalized to one for the 
full distributions, and all slices are plotted in the pair 
center-of-mass (PCMS) frame. Slices are also shown for 
various modifications to the baseline model, in which ex- 
tended time distribution is turned off ("Instant freeze- 
out") and/or oj contribution is removed. The baseline 
model is also shown for the full rj acceptance. The non- 
Gaussian shape in the sidewards direction is due entirely 
to the uj. The same is true for the longitudinal com- 
ponent, where rj acceptance cuts virtually eliminate the 
lifetime effects. The outwards distribution for the uj- 
less, instant freeze-out distribution shows the effect of 
kinematic boost of the PCMS frame, which varies from 
1.7 to 2.8 over the transverse momentum range in this 
bin. The oj and extended lifetime contribute significantly 
above 30 fm, although their combined contributions have 
a less discernible impact. The removal of the rj accep- 
tance cuts has the effect of reducing the impact of these 
contributions along the outwards slice. The angle aver- 
aged distribution shows less sensitivity to the individual 
contributions to the long range source. The angle av- 
eraged sources for the model with and without the ex- 
tended lifetime effects (but with to included) appear to 
be consistent with the PHENIX measurement. 

The dot-dashed line in each plot shows the 1 fm Gaus- 
sian slices for the 3D Bertsch-Pratt fit to the generated 
correlation function (no Coulomb) for the baseline model 
in the PCMS frame, and the bottom panel shows the 
angle-average of this function. Fit results for all model 
variations are tabulated in Table I. The sidewards and 
longitudinal slices show a significant departure from the 
Gaussian shape above 15 fm, yet the outwards distribu- 
tion remains Gaussian in shape out to 40 fm, significantly 
beyond the 30 fm that one would expect from a pure kine- 
matic boost to the PCMS. The sensitivity of the Gaussian 
fits to the outwards is also evident in the radii reported in 
Table I. Beginning with the pure Gaussian input source 
(No UJ, Ar=0) the values of Rside and Riong increase by 
~10% when the lu is included, and Riong increases half 
as much when only the emission tail is added. In con- 
trast, the value of Rout increases by ^60% for either the 
UJ or emission tail alone, and another '^40% when both 
are included. 

In order to investigate these effects in a more complex 
dynamical model with a complete resonance contribu- 
tion with decay channel properties taken from the Par- 
ticle Data Book [39] we show comparable results from 
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FIG. 1: Comparison of the sources sliced in the side, out and 
long directions for each of the cases we consider. The up- 
per panel corresponds to the sidewards direction, the middle 
panel to the outwards direction and the lower panel to the 
longitudinal direction. 



TABLE I: Comparison of the best fit Gaussian parameters in 
each of the cases we consider. We fit directly to the corre- 
lation function in each case. In order to avoid complications 
due to the Coulomb correction, we generated the correlations 
without the Coulomb effect for the purpose of fitting. 



Case 


A 




i?.out 


fl.long 


No uj At = 


0.998 ± 0.006 


3.99 ± 0.02 


5.39 ± 0.03 


4.06 ± 0.02 


At = . 


0.637 ± 0.008 


4.32 ± 0.04 


8.57 ± 0.10 


4.46 ± 0.05 


No u 


0.954 ± 0.008 


3.99 ± 0.03 


8.94 ± 0.06 


4.17 ± 0.03 


Baseline 


0.676 ± 0.010 


4.36 ± 0.05 


12.16± 0.16 


4.45 ± 0.05 


Baseline No Cuts 


0.719 ± 0.007 


4.22 ± 0.03 


6.20 ± 0.05 


4.99 ± 0.04 



the Thcrminator program [40, 41] in Fig. 2. This model 
includes a three-dimensional freeze-out hypersurface sim- 
ilar to the Blast- Wave models [32] including radial flow 
but no emission duration. Particles are created at this hy- 
persurface and propagate freely (without hadronic rescat- 
tering). The model has been shown to reproduce a 
large set of data in the soft sector of the RHIC colli- 
sions [42, 43], and the absolute and relative yields of 
particle species are taken from a thermal model which 
decribes RHIC data well. The distributions shown in 
Fig. 2 are taken from the best-fit parameters [40, 41]: 
T=165.6 MeV, ^lb = 28.5 MeV, r = 8.55 fm/c, p^ax = 
8.92 fm, vt ~ 0.311, using a negative slope in r — t plane. 
The top three panels show 1-fm wide slices along each 
of the Bertsch-Pratt coordinates for the same acceptance 
cuts, \r]\ < 0.35 and 0.2 < kr < 0.36 GeV. For each 
panel we show the full source (all resonances), the full 
source without the oj, and the primordial pion source. 
The full source is normalized to unity, and for the others 
we retain the relative normalization according their con- 
tribution to the two-particle source distributions. For the 
Gaussian slices, we overlay the 1-fm slices from Bowler- 
Sinyukov fits to the full correlation source for a simi- 
lar kx region (average of radii from two bins spanning 
0.15 < kr < 0.35 GeV/c from [41]). To plot the out- 
wards distributions in the PCMS frame, we multiply the 
LCMS Rout fit value by 2.05, corresponding to the kine- 
matic boost for a mean fcy of 0.25 GeV/c. Again we see 
that the outwards slice of the distribution shows the best 
agreement with a Guassian, out to ^38 fm in the PCMS 
(or ~19 fm in the unboosted LCMS), while the slices 
along other Bertsch-Pratt axes deviate from Gaussians 
around 10 fm for longitudinal and 6 fm for sidewards. If 
one projects instead of slices, as was done for Fig. 11 
of [41], the Gaussian shapes extend further in the side- 
wards and longitudinal directions. The angle-averaged 
Thcrminator distribution including all resonances falls 
below the PHENIX image above 20 fm, but accounts for 
~30% of the tail in this region. 

We conclude from this and the previous model com- 
parison that the to is an essential component in the first 
non-Gaussian tail measured at RHIC, but there is suffi- 
cient room for contributions from non-Gaussian lifetime 
effects. The observed sensitivity of the Rout parameter 
to the inclusion of lifetime effects is surprising, and is 
contrary to our initial expectations based on an analysis 
without any experimental acceptance cuts in which the 
Rout parameter showed little change. Only with the ap- 
plication of the PHENIX/STAR 77-cuts did the variation 
appear. This does not rule out the presence of lifetime 
effects, but implies that such affects may already be par- 
tially accounted for in the standard Gaussian fits. Indeed 
the fitted value of Rout for the first model is consistent 
with published PCMS value for Rout for a similar fcy 
range [44]. 

This analysis also underscores the need to perform 
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FIG. 2: Comparison of the Blast- Wave source from the Ther- 
minator model with and without the uj and other resonances. 
We show the angle-averaged sources along with the sources 
sliced in the side, out and long directions. 



source imaging in three dimensions as the most promising 
way to disentangle the to, which has significant contri- 
butions in all directions, from the lifetime effects that 
are restricted to the outwards direction. The source- 
imaging of kaons will provide another, more direct means 
to probe non-Gaussian components of primordial source 
lifetime. Equally important to the task of understanding 
the complete space-time picture will be comprehensive 
set of comparisons to the full source distributions from 
hydrodynamic models with resonances included. 
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